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Abstract 

We develop a noncommutative analogue of the spectral decomposition with the quaside- 
terminant defined by I. Gelfand and V. Retakh. In this theory, by introducing a noncom- 
mutative Lagrange interpolating polynomial and combining a noncommutative Cayley- 
Hamilton's theorem and an identity given by a Vandermonde-like quasideterminant, we 
can systematically calculate a function of a matrix even if it has noncommutative entries. 
As examples, the noncommutative spectral decomposition and the exponential matrices of 
a quaternionic matrix and of a matrix with entries being harmonic oscillators are given. 
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1 Introduction 



The theory of spectral decomposition of a square matrix over a commutative field is well-known 
in linear algebra and is used for calculation of a function of the matrix, especially the exponential 
matrix. However, for a matrix with noncommutative entries, the determinant or the character- 
istic polynomial are not defined because of the ordering problem. Therefore, "eigenvalues" used 
in the spectral decomposition are undefined and we have no systematic method for calculation 
of function of a matrix with noncommutative entries until now. 

Under these circumstances, we studied the exponential of a matrix with entries being har- 
monic osillators for a model in quantum optics and developed "the quantum diagonalization 
method" for a special type of matrices derived from the representation theory [FHKSWJ. More- 
over, we had a chance to encounter with the quasideterminant defined by I. Gelfand and V. 
Retakh. By using the quasideterminant, "noncommutative determinants" such as quaternionic 
determinants [Asj . superdeterminant, quantum determinant, Capelli determainant, etc. are ex- 
pressed in the unified form |GR1] . In the theory of the noncommutative integrable system, 
quasideterminants are very useful to express the solution of the noncomutative integrable equa- 
tions [EGRj, [H], |GNj . |GNUj . Furthermore, various noncommutative analogue of theories using 
determinants are developed, for example, noncommutative analogue of Cramer's formula, the 
Vandermonde determinant, symmetric functions, Pliicker coordinates, and so on. (see, [GR2J, 
[GGRWj . jGKLLRTj and refere nces within). 

In particular, in [ GKLLRTj . they investigated a noncommutative Cayley-Hamilton's theo- 
rem. In their theory, a different characteristic polynomial for each row was introduced and the 
trace or determinant were of the form of diagonal matrices. Moreover, we knew through the 
study of the quantum diagonalization method that eigenvalues should be generalized as "eigen- 
diagonalmatrics" due to the noncommutativity of entries of the matrix. That is why we find that 
a noncommutative Cayley-Hamilton's theorem in [GKLLRTj is suitable to a noncommutative 
analogue of the spectral decomposition. 

In this paper, we define a noncommutative analogue of the Lagrange interpolating poly- 
nomial and develop a noncommutative analogue of the spectral decomposition by using the 
noncommutative Cayley-Hamilton's theorem with the quasideterminant. An identity given by 
a Vandermonde-like quasideterminant plays an essential role. As examples, we explicitly calcu- 
late the noncommutative spectral decomposition and the exponential matrices of a quaternionic 
matrix and of a matrix with entries being harmonic oscillators. 

The contents of this paper are as follows. In section 2, we give a brief review of the spectral 
decomposition in linear algebra. In section 3, we introduce the quasideterminant defined by I. 
Gelfand and V. Retakh and describe some important properties used in our theory. In section 4, 
we review the noncommutative Cayley-Hamilton's theorem in [GKLLRT] shortly. In section 5, 
we develop a noncommutative analogue of the spectral decomposition with the quasideterminant. 
In section 6, we apply our method to a quaternionic matrix and a matrix with entries being 
harmonic oscillators. Section 7 is devoted to disscution. 
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2 Brief Review of the Spectral Decomposition 

Firstly, we give a brief review of the spectral decomposition in linear algebra. 

Let A be a n x n-matrix with commutative entries. For simplicity, we suppose that all the 
eigenvalues Ai, • • • , A n of A are distinct. For j — 1, • • • , n, we set 



n 

l<i<n, iy^j 



(A - \I) 
(A, - A,) 



The polynomial of right hand side is called the Lagrange interpolating polynomial. Then we 
have the spectral decomposition of A; 



A = A1P1 + • • • + \ n P n . 

Moreover, if the Cayley-Hamilton's theorem holdsCthen P±, ■■ 
cesCi.e. 

Pf = P u PiPj^Oii^j), P 1 + --- + P n = I. 
Therefore, we can calculate exp A explicitly; 

exp A = e Xl Pi H h e Xn P n . 

(z - Xi) 



P n are projection matri- 



[Xj Xi) 



U = ,n) 



x 



n—l 



X 



n-l^ 



Remark 1. Lagrange interpolating polynomials fj(z) = J j 

l<i<n, iy^j 

satisfy the following relations; 

(1) x{f 1 (z)+4f 2 (z) + --- + xif n (z)=zi (j = 0,l,...,n-l) 

(2) fi(xj) = 5ij 

We note that from (1), if all distinct, then we have 

\fn(z)J 

3 Quasideterminant 

In this section, we introduce the quasideterminant defined by I. Gelfand and V. Retakh and 
describe some important properties used in our theory. 



-1 


(z n - 


*\ 


) 


[ 1 


I 



3.1 Definition 

Let R be a (not necessary commutative) associative algebra. For a position (i,j) in a square 
matrix A = (a rs )i< riS <„ G M(n,R), let A 1 ^ denote the (n — 1) x (n — l)-matrix obtained from 



A by deleting the z-th row and the j-th column. Let also = (an, 



1 5 



,a in ) and 



c/ = (a lj} 



a V 



Definition 1. We assume that A lj is invertible over R. The (i, j)-quasideterminant of A is 
defined by 



| -A j ij Q>ij 



(1) 
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Example 2. For A 



an &i2 

&21 0,22 



\A\n — an ~~ 021, 



\A\ 2 i = a 2 i - 022^12 an, 



12 — ai2 — ^11^21 «22, 



22 — &22 — 021^11 «12- 



It is sometimes convinient to adopt the following more explicit notation 

\A\ 



11 



an 



021 0,22 



Oil — ftl2«22 °21- 



Remark 3. // the elements aij of the matrix A commute, then 



\A\ij 



_iy+j 



det A 
det A*i ' 



3.2 Inverse Matrix and Quasideteminant 
Proposition 2. If all \A\~j 1 exist, A" 1 is given by 

A' 1 = (\A\ji)l<iJ<n. 



Example 4. For a quaternionic matrix A 



1 % 

j k 



\A\ii 1 = 


(1- 


i ■ 


k- 


■V)- 1 


= (l + 2^r 1 = 


1 

" 2 


\A\2I = 


(J- 


k 


■ % 




= (J + H)" 1 = 


(2J)- 1 = 


\ A \l2 — 


(i- 


1 ■ 


r 




= (^+j'A ; )" 1 = 


m- 1 = 


1 yl 1-1 
l /1 l22 — 


(k- 


j 


■ 1 




= {k-ji)~ 1 = 


(2k)- 1 = 



I 
'2 

i 

2 
k 



Therefore 



A- l = - 
2 



1 -3 
—i —k 



Example 5. We can calculate quasideterminants inductively; 

an - (a 12 a i3 ) 



an 



0-12 Ois 

021 022 023 
031 032 033 



I^ 1 I22 

\ A 1 23 



/1 1 1 1 — 1 
A 1 32 
/1 1 1 1 — 1 
71 1 33 



021 
031 



— an — ^12(022 — a23a 3 3 1 a 3 2) l a 2 i — 012(^32 — 033 a 23 la 22) 1(i 3i 
— <2i3(a23 — a22&Z2 a 3i) 1 021 ~ 013(033 — ^2^22 a ^3,) 1 «31- 
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3.3 Homological Relations 

For A = (a,ij) G M(n,R), n 2 quasideterminants are defined. They are related by the so-called 
homological relations. For example, 



an 


«12 


— — a 22 a 12 1 


an 


Ol2 


021 


a 2 2 


021 


a 2 2 



In general, we have important identities as follows; 
Proposition 3. 1. Row homological relations: 

— I 41 ■ ■ ■ I 4*' I -1 — I 41 ■ I 4^'l _1 

2. Column homological relations: 

— I 4 fc J| -1 ■ I 41 • • — I A i i\~ 1 ■ I 41, • 
I 71 lit I 71 !*? — I 71 Ifct l 71 !^, 



s ^ i 



3.4 The Multiplication of Rows and Columns 

Let B be the matrix obtained from the matrix A by multiplying the i-th row by A G i?Cthen 



\B\ 



* 7 



A|A|jj if fc = z 



if k^i 



Let C be the matrix obtained from the matrix A by multiplying the j-th column by fi G -RCthen 

if / = j ^ 



'^-^ \A\ a if 



Example 6. 



ai2/x 



an 



a 2 i a 2 2/i 



an// 



ai2 



a 2 i£t a 22 
3.5 Sylvester's Identity 



= ai2/i - ana 2 i la 22/U = \A\ 12 fi 
a 12 - a 1 i^(a 2 inY 1 a 2 2 = \A\ 12 



Let A = (a,ij) G M(n, R) be a matrix and A = (oy), i, j = 1, ■ ■ ■ , k, & submatrix of A that is 
invertible over i?. For p, g = fc + 1, • • • , n, set 



A 

dpi ■ ■ ■ dpk 



dig 



"VI 



(3) 



These quasideterminants are defined because matrix A is invertible. 
Consider the (n — k) x (n — k) matrix 

C=(c Pq ), p,q = k + 1,--- ,n. 
The submatrix A is called the pivot for the matrix C. 
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Theorem 4. (Sylvester's identity) For i, j = k + 1, • • • , n, 



\A\ 



1.1 



\C\ij. 



Example 7. 



1 dl2 Ql3 

a 22 a 23 
a 32 a 33 



(i,J = 2,3) 



Applying the theorem H] with the (1, l)-entry 1 as a pivot, we put 

c p<? = n = a pg Cp> 9 = 2, 3) 

and 



1 012 ai 3 
a 22 a23 
a 32 a 33 



|C| 



^22 ^23 
Q>32 «33 



= 2,3). 



4 Noncommutative Version of the Characteristic Polyno- 
mial and the Cay ley- Hamilton's Theorem 

In this section, we review the noncommutative Cayley-Hamilton's theorem in |GKLLRT] shortly. 



We use notations $i(A), Cuu instead of Qi(t), L^\A) in it. For A = ( au ° 12 ) , we denote 

\&21 ^22, 

$i(A), $ 2 (A) as two polynomials given by 



$i(A) = A 2 — (an + ai 2 a 22 a 12 1 )A + (ai 2 a 22 a 12 an — a 12 a 2 i) 
= A 2 - tri(A)X + det i(v4), 



$ 2 (A) = A 2 - (a 22 + a 21 ai 1 a 21 1 )A + (a 21 a n a 21 a 22 - a 21 ai 2 ) 
= A 2 - tr 2 ( A) A + det 2 ( A) . 

Then we can check the noncommutative Cayley-Hamilton's theorem for the generic matrix of 
order 2; 

2 _A ri (A) \ /detail) \ 
^ tr 2 (A)yJ A+ ^ det 2(^)^1 

The general result is as follows. We also give a simple proof. 



Theorem 5. JGKLLRTf For A = (a^-) G M(n,R), we define a "noncommutative characteristic 
polynomial for the i-th row" as follows; 
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*i(A) 



» An) 



"il 



a 



12 



(n-1) (n-1) 
% <2 



a 



(1) a (1) 

(0) JO) 



il 







t2 



= A n — C( ;i ) k \ n 



k=l 



» 



A" 



in , 

,(«-!) \n-l 



'm 

,(0) 



(4) 



where A k = (ajj ). T/ien we /lave a noncommutative version of the Cayley- Hamilton theorem 



I c, 



k=l 



(l)k 



a 



(2)k 



V 



\ 



C(n)fc / 



.4 



n—k 



o. 



(5) 



Proof. For unknown C(i)k (z, fc = 1, • ■ • , n), consider the equation (jSJ). Then the (i, j)-entry of 
©is 



(n-k) 
ij 



(i, j = 1,-- • ,n), 



fc=i 



namely 



(C(i)l' ' " ) C{i)n) 



\ «8> 



•S? / 



a 



(n) 



il i ) "in 



Therefore we obtain C^k by solving the linear equations (171). Moreover, by using 
noncommutative characteristic polynomial for the z'-th row is written as 



(6) 



(7) 



the 



*i(A) 



(n—k) 
il 



(n-1) 

Ai 



(0) 

<1 



En f~i (n—k) 
k=l U (i)fc a in 
>-l) 



A" 



(n-1) 
l il 



A' 



n-1 



,(0) 



o A n - Zfc=i QpfcA 



(n-1) 



n—k 



in— 1 



,(°) 



\n \ v s~i \ n—k 



fe=l 



□ 
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By this proof, we obtain an important corollary. 

Corollary 6. If an identity (G]] holds, the noncommutative characteristic polynomials $«(A) 
defined by (fj]j are equal to X n — J2k=i C{i)k^ n ~ k - Especially, the (usual) Cayley- Hamilton the- 
orem for A (i.e. = Ck for all i) holds, then $j(A) (i = 1, ■ ■ • ,n) coincide with the usual 
characteristic polynomial $(A) of A. 

Moreover, contraposition, we have the following; 

Corollary 7. For a given matrix A, if the noncommutative characteristic polynomials $i(A) are 
different for each i, then no commutative- Cayley-Hamilton-theorem type of identity with respect 
to A exist. 

Example 8. Let A be a matrix A q = (aij) of the generators of the quantum group GL q (n), 
the noncommutative Cayley-Hamilton theorem (the quantum Cayley-Hamilton theorem) holds 
[G KLLRTj . For example, n = 2, by using relations 

a ll a 22 — G 22 a ll = (q 1 ~ q) a 12 a 21i O12O22 = q 1(1 22 a 'Y2i a ll°'21 = <? lfl '21 a 'll> a 12 G! '21 = a 21 ct 12; 

we have 

a2 , 1/2 -1/2 \ fq~ 1/2 \ , , _n , fq- 1 0\ _ 

-Ag - (g 7 a u + Q 1 a 2 2) I Q 1/2 \A q + {a u a 22 - q a 12 a 21 ) I Q J = O. 

However, the noncommutative characteristic polynomials for each row do not coincide each other; 

$i(A) = A 2 - (an + q~ l a 22 )\ + q~ l a X \a 22 - q~ 2 a 12 a 21 , 
$ 2 (A) = A 2 - (q an + a 22 ) A + q a u a 22 - a 12 a 2 i. 

Therefore, there is no identity for A of commutative-Cayley-Hamilton-theorem type. 

5 Noncommutative Spectral Decomposition 

In this section, we develop a noncommutative analogue of the spectral decomposition with the 
quasideterminant. First, we review the Vandermonde quasideterminant and define a noncommu- 
tative analogue of the Lagrange interpolating polynomial. Next, we present the main theorem 
and our method of a noncommutative spectral decomposition. We also give a proof of the 
theorem by using properties of the quasideterminant prepared in section [HJ 

5.1 Vandermonde Quasideterminant 

First, for xi, x 2 , ■ ■ ■ , x k G R , the Vandermonde quasideterminant ( |GKLLRTj . [GGRW]) is 
defined by 









x k 


V(x ir - 


■ ,x k ) = 


X\ 

1 


Xk 
1 
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Example 9. 

V(xi,X2,Z) 



z 



™2 ™ 2 

1 2 



Xl X2 Z 
111 



2 _ / 2 2W ^1 ^2 

- l x i x 2) l( X2 - Xl )-i (l-x^)- 1 



= -2 2 + (-xi - (x 2 - xi)x 2 (x 2 - xi) 1 (x 2 - xi)x 2 (x 2 — Xi) 1 x 1 ) L 

= ^ 2 + + 2/2) Kj 

= z 2 - (j/i + y 2 )z + y 2 yi 

where we put y\ = Xi, 1/2 = (x 2 — Xi)x 2 (x2 — xi) -1 . This is the noncommutative version of 
the relationship between solutions and coefficients for a (left) algebraic equation of degree 2 
[GtCtRW] . 



Remark 10. If z = A = (a*,,-) G M(2,R) and Xj 



X(i)j 



x 



(2)j 



(j = 1,2), j/j- (j = 1,2) are 



a/50 diagonal matrices. Moreover, comparing the equation V(x±, x 2 , A) = A 2 — (yi+y 2 )A+y 2 yi 
O /ne noncommutative Cayley-Hamilton's theorem 



A - 



C ( i)i 



C(2)l 



a 



(2)2 



«/ 2/i + 2/2 



1 " „ ) and 2/22/1 = — ( „ ) , fry i/ie relationship between solutions 

G (2)i/ V G (2)2/ 

and coefficients again, X(iU,X(j) 2 are the solutions of the noncommutative characteristic equation 

of A for the i-th row. 

For a given z = A = (a^) G M(n, R) and the equation V(xi, • • • , x n , A) = O, diagonal com- 
ponents of diagonal matrices Xj are the solutions of the noncommutative characteristic equations 
of A in the same way. 

5.2 Noncommutative Lagrange Interpolating Polynomial 

(n-l ™n-l N 
1 .-. 1 

exist. Then we define polynomials fc(z) (i = 1, • • • , n) with respective to z G -R as follows; 
Definition 8. 



(h 




-( 


\fn 







_n-l 



X 



n-l N 



V 1 / 



We call them noncommutative Lagrange interpolating polynomials. 
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Example 11. For n = 2, 



| Xi 


x 2 


-1 

z + 


X\ 


x 2 


1 


1 


1 


1 



(xi - x 2 y 1 z + (1 - x 2 l X\)~ l 

(x 1 - x 2 y 1 z + (x 2 - x 1 y 1 x 2 

(xx - x 2 )~ l (z - x 2 ), 



Hz) 



X\ x 2 

1 1 



z + 



Xi x 2 

l VT 



x 2 — X\) 1 z + (1 — x 1 l x 2 ^ 1 
1 



= [X 2 -Xi) [Z-Xi). 

By the definition above, we obtain the following theorem. 

Theorem 9. For xi, ■ ■ ■ , x n , z G R, we have 

(1) x\h(z) + xif 2 (z) + ■ ■ ■ + xlfn(z) =2? (j = 0,l,...,n-l), 

(2) friz,) = Sij. 

5.3 Our Method of Noncommutative Spectral Decomposition 

Theorem 10 (Main theorem). If given z, xi, ■ ■ ■ ,x n G R satisfy the equation V(x\, • • • , x n , z) 
0, then we have the following identities 



Vrr, 



X 



X 



1 

n-1 



rf.n—1 ~n—l 



1 



(m = 0, • • • , n, n + 1, 



Rewriting (jHJ), by the definition of noncommutative Lagrange interpolating polynomials 



„n-l 



X i 



.r 



m \ 
n ) 



1 




-1 

1 


(z n - 


!\ 








{ 1 


/ 



= *r/i(*) +•••+</»(*). 

then we have the noncommutative spectral decomposition of 2; 

z m = + • • • + x™f n (z) (m = 0, 1, • • • )• 

In particularCif 2 is a matrix A = (a^) G M(n, i?), putting xi, • • • ,i n as unknown diagonal 
matrices and solve the equation 

V(xx, ■ ■ ■ ,x n , A) = O. 

By the remark [TU1 this equation is nothing but the noncommutative Cayley-Hamilton's theorem 
and the diagonal components of diagonal matrices Xj are the solutions of the noncommuta- 
tive characteristic equations of A. Therefore, by using the solutions of them, we obtain the 
noncommutative spectral decomposition of A 

A m = xTh(A) + ■■■ + x™f n (A) (m = 0, 1, • • • )• 
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5.4 A Proof of Main Theorem ITTJ 



Proof. In case oi m — 0,1, ■ • ■ ,n — 1, the identity (El) is trivial. If m 
V(xi, ■ ■ ■ , x n , z) = 0. In the following, we suppose m — n + 1, n + 2, • ■ • 
Consider a matrix ^4 and the submatrix defined by 



n, (E]) is nothing but 



/ 



.4 



.71 



v 



^ 



.r 



n-l 



7 



For p = 1,2, g = n + 1, n + 2, we put a matrix C = (c P9 ) which entries are quasideterminants 
with A as a pivot like ([3]) (note that quasideterminants are unchanged under permutations of 
rows or columns) and we remark 

... rr-n 



Cl,ra+2 



^ rr, 



v n— 1 



Then by the Sylvester's identity (theorem H]), we have 



\A\l,n 



+2 



\C\l,n 



+2 



Cl,n+1 Ci )Tt+ 2 



C2,n+1 C2 )Tt +2 
1 



Cl,n+2 — Ci jn _|_iC 2jn+1 C2 in +2 
V m ~ c l,n+l c 2,n+l^- 
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On the other hand, since 





x m 


rf-m 
x n 





z m 




x l 


■^n 


2 m 







X-^ 


x n 





z n 






•^n 







|^|l,n+2 = 


X 


rpTl— 1 





z n-i 






~n— 1 


z n-l 







1 


1 


1 


1 




1 


1 


1 


1 



X 



X 



rn 
1 

X-y 
.n-1 
1 



m— 1 



.r 



rn 
n 
n 
n 
n-1 
n 



Z 

„n-l 



1 

n— 1 
1 

„n-2 



a; 
.x i 



m— 1 



by the Sylvester's identity with 
(rt + 2, n + 2)-entry 1 as a pivot 



n 

n-1 

n 

n-1 
n 



y n—l 
v n-2 



(by the property ([2])) 



1 

we obtain the important identities 

^m— \Z — V m c l,n+l c 2,n+l^n- 

Therefore V n = implies = (m = n + l,n + 2,---)by the mathematical induction. 
Remark 12. In particular, for n = 2, 



V m = z m - (x? x%) (J ^ 
T/ien i/ie identity (TJ|) is 



z m -{x™f x {z)+x™f 2 {z)). 



„m— 1 „,m— 1 N 



(9) 
□ 



(Xn 



x?- x ){x 2 - x 1 y 1 {z 2 - (xiMz) + xif 2 {z))} (m = 2, 3, • • • ). (10) 



6 Examples of Noncommutative Spectral Decomposition 
and the Exponential Matrices 

In this section, we apply our method to a quaternionic matrix and a matrix with entries being 
harmonic oscillators. As a result, we obtain the noncommutative spectral decomposition and 
the exponential matrices of them explicitly. 
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6.1 Quaternionic Matrices 

As a quaternionic matrix, we consider an element A of Lie algebra sp(2)\ 

A={\ 3 
v ~h 

We apply our method to A and calculate the spectral decomposition and the exponential matrix 
exptA explicitly. First, from 

A 2 _( -2 2k 
~ {-2k -2 / 

noncommutative characteristic equations for each row are 

-2 2k fA^ 



then A = 0, 2i, and 



$i(A) 



$ 2 (A) 



3 A 
1 



A - 2i\ = 0, 



-2k ! 





2 [T[ 
? A 
1 1 



A 2 + 2iA = 0, 



then A = 0, -21 

Next, in the noncommutative Lagrange interpolating polynomials 

= (x x - X 2 y 1 {z - X 2 ), f2(z) = (X 2 - XiY 1 (z - Xi), 



we put z = A, 
Pi 



2i 



fi(A) 



MA) 



-2% 
'2% 



x 2 











then 



-2% 
2i 



-i 



A 



1 fl 



-2i 



2 \k 



A 



2i 



-2% 



1 k 
-k 1 



We can check Pf = Pi, PiPj = (i ^ j) easily and we obtain 



-,2it 



exptA = (exp txi)Pi + (exp tx 2 )P 2 

1 / 1 —k 

2 \k 1 

-e 2it k + k 



1 



e 

2it 
2it, 



-lit 



+ 



1/1 fc 



1/ 2 V-lfe 1 



+ 1 



2 \e-^fc - k 



-2it 



+ 1 



Remark 13. If we put 



e cos t e sin t j 
e~ rf sin t j e~ lt cos t 

2% 



e Sp{2). 



X\ 







x 2 







2/ 



we have [A, xi] ^ O, [A,x 2 ] ^ O and fi(A), f 2 (A) are not projection matrices. Nevertheless, 
by the theoremUM we have A m = x™fi(A) + x^f^A) and we can calculate expL4 explicitly! 
This result is derived from the fact that the theoremlWiis not depend on the ordering of solutions 
for noncommutative characteristic equations for each row. 
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6.2 A Matrix with Entries being Harmonic Oscillators 

Let a, a) be the generator of the harmonic oscillatorD The relation is [a, a*] — 1. We also denote 
N as the number operator TV = a) a. 

/0 a 0\ 

We consider a matrix A = \/2 I a j . This matrix is related to a Hamiltonian of a 

\0 at 0/ 

model in quantum optics [FHKSWJ. So, it is important to calculate the exponential of A as the 
time-evolution operator of the Hamiltonian. 
From 



A 








2N + 1 , A 3 = 2V2 (2N + 1) 











AT, 







(2iV + 3)a 


(2N- l)at 





(27V+l)a 




the noncommutative characteristic equations for each row are 
$i(A) 







2(JV+1) 

1 



2x/2(2iV + 3)a [A^ 

2a 2 A 2 " 
v^a A 

1 



A 3 - 2(2N + 3)A = 0, 



then A = ±^2(2^ + 3), 0, and 
$ 3 (A) -- 



2 v ^(2iV-l)a t [F 

2(at) 2 2NX 1 

v^at A 

11 



A 3 - 2(2iV- 1)A = 0, 



then A = ±y/2(2N-l), 0. 

Remark 14. For the second row, the quasideterminant 

2V2(2N+l)a) 2^(2^ + 1)0^ 



$ 2 (A) 






2(2JV+ 1) 





A 2 







V2a 


A 





1 





1 



2(2JV+1) 

is not defined because "rank" of the matrix ( y/2a) \[2a is 2 (on "rank of A E 











M(n, R) ", see IGGRWf ). Then, we put 

A 3 + UA 2 + VA + W = 0, U, V, W are diagonal matrices 
and simplify them, then the second row is 

A 3 + u\ 2 - 2(2N + 1)A - 2u(2N + 1) = (for arbitrary u). 



Therefore, if we put u = 0, we Main A = ±y/2(2N+l), 0. 
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Next, we calculate the noncommutative Lagrange interpolating polynomials 



h{z) 



X 



tjC o t]C q 



Xi x 2 x 3 
1 1 1 



z 2 + 







-1 












x 2 


x 3 


z + 


Xi 


1 


1 




1 



= Xi 



xl(x 2 -x 3 ) 1 xi - x\(x 3 - x 2 ) l xi 



-x\(x 3 -x 2 ) 1 x 3 - x\(x 2 - x 3 ) l x 2 ) l z 2 



+ (xi - x 2 (x 2 - x 3 ) x 1 -x 3 (x 3 -x 2 



-x 2 (x 3 - x 2 ) 



x 3 (x 2 - x 3 ) x 2 ) z 



+ (1 - (x 2 - x 3 x 2 ) x l - (x 3 - x 2 x 3 ) x l 

-(X 2 - X 3 1 X 2 2 )~ l Xi - (x 3 - X^xD^Xi)' 1 . 

In particular, in the case of x\ = x, x 2 = 0, x 3 = —x, 

fi(z) = (x 2 -(-xfi-x)-^)-^ 2 + (x-(-x)(-xy 2 x 2 )~ 1 z + 
= (2x 2 y 1 z 2 + (2x)- 1 z, 

where the last term of fi (z) is calculated by using the homological relation as follows; 



2 2 2 

ry* /y* ry* 

x l x 2 x 3 

Xl x 2 x 3 



1 



2 2 2 
1 2 3 



Xl I x 2 x 3 

1 1 1 



-1 



2 

x l 2 


2 

x 2 3 




x\ 


x\ 


X 2 


x 3 




1 


1 



*^ 1 



Xi 



x 2 x 3 

1 1 1 



(X 2 - X 3 1 X2)(1 - X 3 2 xl) 1 

(2x) _1 -0-1 = (xi — > x, x 2 — > 0, x 3 — > — x). 



Then we put 



A, 



x 



A(iV - 1) 



A(iV-2), 



A(iV) = v/2(2iV + 3), 
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we have 



Pi = fi(A) 

= (2x 2 Y l A 2 + (2X)- 1 A 
f(2X(N))- 2 

i-2M.\ 1 )) 2 I -2 

\ (2X(N-2))- 2 

\2XiN))- 1 \ /O a 

- ? (2A(iV - l))- 1 • V2 ^ a 

(2A(7V-2))-y \0 at 



iV + 1 a 2 
2iV + 1 
TV, 



(at) 2 



/ 1 

' 2(2iV+3) 



2(2iV+l) 



2(2iV-l)/ 



/ 1 

' 2V27V+3 



+ 



^ /iV + 1 a 2 \ 
2N + 1 



i 



2V27V+1 



2V27T 



) 2 









/o 


a 


0\ 




at 





a 




VO 


at 


0/ 



JV+l 
2(27V+3) 

^^at 



2 v / 2TV+3~ a 2(2iV+3) 
1 1 



a 2 \ 



2 % /27V+T 2 2V27V+1 

^2(2AT-1) ( flt ) 2 2^2^-1 flt 2(27V-1) / 



In the same manner, if we put z = A, x\ — x, x 2 = 0, x 3 = —x in f2(z), fs(z), then we have 
P2 = f 2 (A) 



(-x 2 )- l A 2 + h 



/ l 

' 2(27V+3) 



2(27V+1) 



\ 



2(2iV-l)/ 





(N + l 





a 2 \ 




n 


\ 


2 





IN + 1 





+ 




i 




I (« f ) 2 





NJ 




\ 


V 



7V+2 



2N+3 










2N 



2N+3 

h-Mf o 1 



1 a 2N 



2JV-1 
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P 3 = f 3 (A) = (2x 2 )- 1 A 2 -(2x)- 1 A 



( 



l 



2(2iV+3) 



2(27V+l) 



1 



2(2JV-1)/ 



'N + 1 a 2 \ 
2iV + 1 
JV/ 



a 



t^2 



2V27V+3 



2V277+T 



/ 



jV+1 
2(2iV+3) 
1 



1 



-a 



2V2W+1 

\2(2JV-1) ^ U I 



2V21V+3 

2 
1 



2y/2N-l/ 
\ 







a 


0\ 




at 





a 


1/ 


\o 


at 


0/ 



2(2iV+3) 
1 



2V27V-1 



2 v / 2lV+T 
U 2(2AT-1) / 



We can check Pf = P iy P%Pj = (i 7^ j) and for a constant g, 

exp(-itgA) 
= exp(—itgx)Pi + (exp 0)p2 + exp(ztc/:r)P3 



/ A r +2+(A f +l) cos(tg\(N)) 

' 2N+3 

-i^jsin(^A(JV-l))ot 



i^— sin(t</A(iV))a 
cos(tg\(N- 1)) 



-1 + cos(tgA(iV)))a 2 \ 



2N+3 

- i ^+T sin (^A(iV-l))a 

N-l+N cos(tg\(N-2)) i 
2N-1 / 



VWr T (- 1 + cos(t^A(iV - 2)))(at) 2 -^ 73 ^ T sin(t 5 A(iV - 2))at 
Remark 15. For i/iis A, by using "the quantum diagonalization method" r L FHKSW], 



ft 



V 



\ 



A 



y/N(N-l)/ 



ft 



\ 



\ 



VN 







y/N(N-l) y ) 



VNTT 

V2 1 v^vTI v^VT2 



VN + 2 







Since the matrix on the right-hand side has only commutative entries, we calculate the charac- 
teristic equation as usual, then 

A 3 -2(2iV + 3)A = 0, A = 0, ±^2(2^ + 3). 

We remark that the result of exp(-itgA) in [FHKSWJ and with our noncommutative spectral 
decomposition coinside. 



7 Discussion 

In this paper, we developed a noncommutative version of the spectral decomposition with the 
quasideterminant and calculated some interesting examples. In particular, we defined a non- 
commutative analogue of the Lagrange interpolating polynomials and applied to the systematic 
method for constructing projection matrices with noncommutative entries. 

Our method is very powerful to calculate a function of a matrix with noncommutative entries 
and is expected to apply for the theory of noncommutative geometry, quantum physics, and so 
on. A study of other applications with our theory is in progress. 
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